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Abstract: We describe the application of a novel approach for the reduction of scatter- 
ing amplitudes, based on multivariate polynomial division, which we have recently pre- 
sented. This technique yields the complete integrand decomposition for arbitrary ampli- 
tudes, regardless of the number of loops. It allows for the determination of the residue 
at any multi-particle cut, whose knowledge is a mandatory prerequisite for applying the 
integrand-reduction procedure. By using the division modulo Grobner basis, we can derive 
a simple integrand recurrence relation that generates the multi-particle pole decomposition 
for arbitrary multi-loop amplitudes. We apply the new reduction algorithm to the two- 
loop five-point planar and non-planar scattering amplitudes in J\f = 4 SYM and A^ = 8 
SUGRA in four dimensions, whose numerator functions contain up to rank- two terms in 
the integration-momenta. We determine all polynomial residues parametrizing the cuts of 
the corresponding topologies and sub-topologies. Moreover, we obtain the integral basis 
for the amplitude decomposition, defined by the polynomial form of the residues. Our 
approach is well suited for a semi- numerical implementation, and its general mathemati- 
cal properties provide an effective algorithm for the generalization of integrand reduction 
method to all orders in perturbation theory. 
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1. Introduction 

The unitarity of the (S-matrix encodes the most profound property of a quantum system, 
namely the probability conservation. The optical theorem, that relates the difference be- 
tween the transition amplitude and its complex conjugated one to their product, is the 
direct consequence of unitarity. Hence, at a given order in perturbation theory, it connects 
the discontinuity of the amplitude across a given branch cut to the sum of all the Feynman 
diagrams sharing that specific cut, which factorize into two lower-order amplitudes. By 
elaborating on the role of the optical theorem, and introducing the concept of generalised 
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cuts [1—4], unitarity has been inspiring a novel organization of the perturbative calcu- 
lus, where Feynman diagrams are grouped according to their multi-particle factorization 
channels. 

Scattering amplitudes in quantum field theories are analytic functions of the momenta 
of the interacting particles, hence they are determined by their singularities. The singu- 
larity structure is retrieved when virtual particles go on-shell, under the effect of complex 
deformations of the kinematic, as needed for solving multiple on-shell conditions simulta- 
neously. 

The investigation of the mathematical properties of the residues at the singularities 
led to the discovery of new relations involving scattering amplitudes, such as the BCFW 
recurrence relation [4], its link to the leading singularity of one-loop amplitudes [2], and 
the one- loop OPP integrand-decomposition formula [5]. 

Automating the evaluation of one-loop multi-particle amplitudes, for an accurate de- 
scription of scattering processes that were considered prohibitive, has become feasible. 
Motivated by the challenging experimental programme of the LHC, where the ubiquity of 
QCD manifests itself through the production of multi-jet events, several codes have been 
developed with the goal of reaching the next-to-leading order level of accuracy for the 
cross-sections [6-13]. 

On the more mathematical side, it became clear that within the on-shell and unitarity- 
based methods, the theory of multivariate complex fuctions could play an important role 
in order to compute the generalized cuts efficiently. The holomorphic anomaly [14, 15] 
and the spinor integration [16,17], as well as, Cauchy's residue theorem [2,4], Laurent 
series expansion [18-20], Stokes' Theorem [21,22], and Global residue theorem [23] have 
been employed for carrying out the integration of the phase-space integrals, left over after 
applying the on-shell cut-conditions to the loop integrals. 

Progress on the unitarity-based methods and the vivid research activity spinned off 
has been recently reviewed in [24] and [25-31]. 

At two loops, generalized unitarity techniques have been introduced for supersym- 
metric amplitudes [32] and later for QCD amplitudes [33]. The multiple cuts of two- loop 
amplitudes were proposed to extend the simplicity of the one-loop quadruple-cut [2] to the 
leading singularity techniques [34,35] and to the maximal cuts method [36]. The maximal- 
unitarity approach developed by Kosower, Larsen, Caron-Huot, and Johansson [37-40] has 
refined this techniques by a systematic application of the global residue theorem. 

The singularity structure of multi-loop scattering amplitudes can be also exposed in 
the integrand, before integrating over the loop momenta. The integrand reduction methods 
use the singularity structure of the integrands to decompose the (integrated) amplitudes 
in terms of Master Integrals (Mi's). The multi-particle pole expansion of the integrand is 
equivalent to the decomposition of the numerator in terms of products of denominators, 
multiplied by polynomials. These latter correspond to the residues at the multiple-cuts. 
In general, the coefficients of the Mi's are a subset of the coefficients appearing in the 
polynomial residues. Therefore the complete determination of the residues leads to the 
complete decomposition of the amplitudes in terms of Mi's. 
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The parametric form of the polynomial residues is process-independent and it can be 
determined a priori, from the topology of the corresponding multiple cut, namely from the 
diagram identified by the denominators that go simultanously on-shell. The actual value of 
the coefficients is clearly process-dependent, and its determination is indeed the goal of the 
integrand reduction. Integrand reduction methods determine the (unknown) coefficients 
by polynomial fitting, through the evaluation of the (known) integrand at values of the 
loop-momenta fulfilling the cut conditions. 

The integrand reduction methods have been originally developed at one loop [5]. Ex- 
tensions beyond one- loop were proposed in [41,42]. A key point of the higher- loop extension 
is the proper parametrization of the residues of the multi-particle poles. Each residue is 
a multivariate polynomial in the irreducible scalar products (ISP's) among the loop mo- 
menta and either external momenta or polarization vectors constructed out of them. ISP's 
cannot be expressed in terms of denominators, thus any monomial formed by ISP's is the 
numerator of a potential MI which may appear in the final result. 

Both the numerator and the denominators of any integrand are multivariate polynomi- 
als in the components of the loop variables. As recently shown in [43,44], the decomposition 
of the integrand can be obtained by performing the multi-variate polynomial division be- 
tween the numerator and the denominators, using basic principles of algebraic geometry, 
like the division modulo Grobner basis. Moreover, the multi-variate polynomial divisions 
gives a systematic classification of the polynomial structures of the residues, leading to 
both the identification of the Mi's and the determination of their coefficients. 

In [41] it was observed that the set of independent integrals which emerge from 
the integrand reduction algorithms is not minimal. Integration-by-parts [45], Lorentz- 
invariance [46], and Gram-determinant [47] identities may constitute additional, indepen- 
dent relations which can further reduce the number of Mi's that have to be actually eval- 
uated, after the reduction stage. Badger, Frellersvig and Zhang have explicitly shown 
that the number of independent 10-denominator independent integrals identified through 
the integrand decomposition of the three-loop four-point ladder amplitudes is significantly 
reduced by using integration- by-parts identities [48]. In the case of one- and two- loop 
amplitudes, an alternative technique for counting the numbers of tensor structures and of 
the independent coefficients has been presented by Kleiss, Malamos, Papadopoulos and 
Verheynen in [49]. During the completion of this work, Feng and Huang [50] have shown 
that, by using multivariate polynomial division [43,44], a systematic classification of a four 
dimensional integral basis for two-loop integrands is doable. 

In [44] , we have set the mathematical framework for the multi-loop integrand reduction 
algorithm. We have shown that the residues are uniquely determined by the denominators 
involved in the corresponding multiple cut. We have derived a simple integrand recur- 
rence relation generating the multi-particle pole decomposition. The algorithm is valid for 
arbitrary amplitudes, irrespective of the number of loops, the particle content (massless 
or massive), and of the diagram topology (planar or non-planar). Interestingly, at one- 
loop our algorithm allows for a simple derivation of the one-loop OPP reduction formula. 
The spurious terms, when present, naturally arise from the structure of the denominators 
entering the generalized cuts. 
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In the same work [44], we gave the proof of the maximum-cut theorem. The theorem 
deals with cuts where the number of on-shell conditions is equal to the number of inte- 
gration variables and therefore the loop momenta are completely localised. The theorem 
ensures that the number of independent solutions of the maximum-cut is equal to the num- 
ber of coefficients parametrizing the corresponding residue. The maximum cut theorem 
generalizes at any loop the simplicity of the one- loop quadruple-cut [2,5], where the two 
coefficients parametrizing the residue are determined by the two solutions of the cut. 

In this paper, we apply our algorithm to the two-loop five-point planar and non-planar 
amplitudes in M = 4 super Yang-Mills (SYM) and M = 8 Supergravity (SUGRA) in four 
dimensions [51,52], whose numerators contain up to rank- two terms in each integration- 
momenta. In particular, we derive the generic polynomial residues which are required by 
the reduction procedure. Later, we show that the integrand reduction can be performed 
both semi-numerically, by polynomial fitting, and analytically. The latter computation 
has been performed generalizing the method of integrand reduction through Laurent ex- 
pansion [53], which has been recently introduced to improve the integrand reduction of 
one-loop amplitudes. 

All the numerical and analytic computations presented in this paper have been per- 
formed using C++, FORM [54] and the mathematica package S@M [55]. 



2. Integrand reduction 

In this Section we describe the general strategy for the reduction of scattering amplitudes 
at the integrand level, following [41,44]. In dimensional regularization an £-\oop amplitude 
with n denominators is written as 

An = J d d qi ... j d d q e Ti v -i n (qi, ■ ■ ■ , q~e) 

= fd d q!... fd% -A/ii-Ugi. ••-.&) 

J 1 "J Ai (qi, ■■■,&)■■■ A„(?i, ■■■,qe) ' 

A = ^2 ai, a q a + p}j 



mf (2.1) 



where qi,...,qg are integration momenta and aj j0 € {0, ±1}. Objects living in d = 4 — 2e 
are denoted by a bar. We use the notation qa = qa + \ a , where qH is the four-dimensional 
part of q a , while X qa is its (— 2e)-dimensional part [56]. In the following we will limit 
ourselves to the four-dimensional case. Extensions to higher-dimensional cases, according 
to the chosen dimensional regularization scheme, can be treated analogously. 

The integrand reduction methods [5,41,42,53,57-64] trade the decomposition of the 
loop integrals in terms of Master Integrals (Mi's) with the algebraic problem of building a 
general relation, at the integrand level, for the numerator function of the amplitudes. In 
this paper we use the method introduced in [44]. The algorithm relies solely on general 
properties of the loop integrand, i.e. on the maximum power of the loop momenta present 
in the numerator, and on the quadratic form of Feynman propagators. The residue of 
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each multi-particle pole is determined by the on-shell conditions corresponding to the 
simultaneous vanishing of the denominators it is sitting on. In particular, we obtain the 
multipole decomposition of the integrand of Eq. (2.1) using an integrand recurrence relation 
based on multivariate polynomial division together with a criterion for the reducibility of 
the integrands. In the following subsections we will briefly review the two ingredients of 
the method. 

2.1 Integrand recurrence relation 

The 4-dimensional version of the integrand of Eq. ( |2.1| ) is 

2 -A/ii-i^gi, ■■■,?<) ,„ „x 

iv " ln D il (q 1 ,...,q e )---D tn (q 1 ,...,q e ) ' 1 ' > 

The numerator J\fi 1 ...i n and any of the denominators A are polynomial in the components 
of the loop momenta, say z = {z\, . . . z±e), i.e. 

■j _ ^il-»n( Z ) /r, o\ 

Ai (*)••■ A„(z) 

We construct the ideal generated by the n denominators 

Jn-i n = (Ax,-- - ,AJ = |X>k(z)A>) 

where P[z] is the set of polynomials in z. The common zeros of the elements of Sii—i n 
are exactly the common zeros of the denominators. We chose a monomial order and we 
construct a Grobner basis generating the ideal Ji x ...i n 

g il ... in = {g 1 (z),...,g m (z)} . (2.4) 

The n-ple cut-conditions D. L1 = . . . = Di n = 0, are equivalent to g\ = . . . = g m = 0. The 
multivariate division of Mi v ..i n modulo Qi x ---i n leads to 

M 1 ... in (z)=r il ... i „+A il ... in (z) , (2.5) 

where r^...^ = Y^i=\ Qi( z )9i{ z ) i s a compact notation for the sum of the products of the 
quotients Qi and the divisors gi. The polynomial A^...^ is the remainder of the division. 
Since Q% x —i n is a Grobner basis, the remainder is uniquely determined once the monomial 
order is fixed. The term T^...^ belongs to the ideal J7ix-i n , thus it can be expressed in 
terms of denominators, as 

n 

IV-in = ^M r .. iK _ lWl ... in (z)A K (z) . (2.6) 

K=l 

The explicit form of Mi v ..i K _ 1 i K+1 ...i n can be found by expressing the elements of the Grobner 
basis in terms of the denominators. Using Eqs. ( |2,5| ) and (|2.6| ), we cast the numerator in 
the suggestive form 

n 

Mi-i» = X^-u-iWi-inC^A^z) + A n ... in (z) . (2.7) 

K=l 
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Plugging Eq. ( |2.7D in Eq. (2.3), we get a non-homogeneous recurrence relation for the 
n-denominator integrand, 



A, 



■Eii—in ~ ^~]-^h—i K -ii K +iin ^ n. . . . n. ' $ 



According to Eq. Zi 1 ---i n is expressed in terms of (n — l)-denominator integrands, 

q- _ ■A/"i 1 ...t K _ 1 i K + 1 j n 

- ^ _ _ _ d . k _ iD . . . D . n ■ W 



The non-homogeneous term contains the remainder of the division fl2.5|) . By construction, 
it contains only irreducible monomials with respect to Gi v -i n , and it is identified with the 
residue of the cut {i\...i n ). 

The integrands Ii 1 ---i K _ 1 i K+1 -i n can be decomposed repeating the procedure described in 
Eqs. ( p.3[ )-(^). In this case the polynomial division of J\fi 1 —i K _ 1 i K+1 ...i n has to be performed 
modulo the Grobner basis of the ideal Ji 1 —i K _ 1 i K+1 —i n , generated by the corresponding 
(n — 1) denominators. The complete multi-pole decomposition of the integrand 2^...^ is 
obtained by successive iterations of Eqs. (|2 . 3|) - (|2 . 5[) . 

2.2 Reducibility criterion 

An integrand 2j r ..j n is said to be reducible if it can be written in terms of lower-point 
integrands, i.e. when the numerator can be written as a linear combination of denominators. 



Eqs. ( |2.5|) and (2J3), allow to characterize the reducibility of the integrands: 



Proposition 2.1 The integrand 1^...^ is reducible iff the remainder of the division modulo 
a Grobner basis vanishes, i.e. iff A/i r ..j n S Ji\-i n - 



A direct consequence of the Proposition 2.1 



is 



Proposition 2.2 An integrand Ii 1 ...i n is reducible if the cut {i\ ■ ■ ■ i n ) leads to a system of 
equations with no solution. 

Indeed if the system of equations (z) = • • • = Di n (z) = has no solution, the weak Null- 
stellensatz theorem ensures that 1 E Ji x ...i n , i.e. Ji x ---i n = -P[z]- Therefore any polynomial 
in z is in the ideal. Any numerator function Afi l ...i n is polynomial in the integration mo- 
menta, thus Mi v ..i n E Ji± - i n and it can be expressed as a combination of the denominators 
D h (z), . . . , D in (z) [44, 49]. In this case Eq. ([T|) becomes 

n 

■^il—in = y ; -^H-"i«-lWl»-ra ' (2-10) 

K=l 

The reducibility criterion and the recurrence relation (|2.8| ) are the two mathematical 
properties underlying the integrand decomposition of scattering amplitudes, at any order 
in perturbation theory. If the n denominators cannot vanish simultaneously, the corre- 
sponding integral is reducible, namely it can be written in terms of integrands with [n — 1) 
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denominators. If the n-ple cut leads to a consistent system of equations, we extract the 
polynomial form of the residue as the remainder of the division of the numerator modulo 
the Grobner basis associated to the n— pie cut. The quotients of the polynomial division 
generate integrands with (n — 1) denominators which should undergo the same decompo- 
sition. The algorithm will stop when all cuts are exhausted, and no denominator is left. 
Upon integration, the non-vanishing terms present in each residue may give rise to master 
integrals. 

Each residue A^...^ belongs to a vector subspace Qi r ..i„[z] of P[z]. Its dimension is 
independent of the choice of the basis. The residue can be expressed in terms of the ISP's 
writing the components z in terms of scalar products. A suitable choice of the bases of the 
loop momenta allow to write the ISP's as multivariate monomials in z generating Qi v ..i n [z]. 
When the number of external legs of the cut diagram is less than five, then the ISP's may 
involve spurious terms. As in the one-loop case [5,41,42,53,57-63], they originate from 
the components of the loop momenta belonging to the orthogonal space, i.e. the space 
orthogonal to the one spanned by the independent external momenta of the cut diagram. 

2.3 Maximum-cut Theorem 

Following [44], we define Maximum cut an n-ple cut Di 1 (z) = ■■■ = Z)j n (z) = fully 
constraining all the components z of the loop momenta. Examples of maximum cuts for 
the one-loop case are the 4-ple (5-ple) cut in 4 [d = 4 — 2e) dimensions. We assume 
that, at non-expceptional phase-space points, a maximum cut has a finite number n s of 
solutions, each with multiplicity one. Under these assumptions one can prove the following 
theorem [44] 

Theorem 2.1 (Maximum cut) The residue at the maximum-cut is a polynomial para- 
matrised by n s coefficients, which admits a univariate representation of degree (n s — 1). 



The maximum-cut theorem guarantees that the maximum number of terms needed to 
parametrize the residue of the maximum cut is exactly equal to n s . Therefore it guarantees 
the full reconstruction of the residue by sampling the integrand on the n s solutions of the 



maximum cut. The Theorem |2.1| generalizes at any loop the simplicity of the one-loop 
maximum cuts [2,5]. Indeed, in d = 4 — 2e dimensions, the residue of the quintuple cut is 
parametrized by 1 coefficient, and can be reconstructed by sampling on the single solution 
of the cut itself. Similarly the two coefficients of the residue of the quadruple-cut in 4 
dimensions can be determined by sampling the integrand on the two solutions of the cut. 

2.4 Two-loop integrand reduction 

In 4 dimensions the generic 2-loop n-denominator amplitude is 

N\... n {q,k) 



A n = j d A q J d 4 k l n (q,k) = J d 4 q J d A k 



Di D 2 ■ ■ ■ D n ' 

A = (ai,i<? + a 2 ,ik + pi) 2 - m 2 (2.11) 
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Every integrand with more than eight denominators -D, leads to a system of equations with 
no solution for its cut 1 . Proposition 2.2 implies that such an integrand is reducible and 
can therefore be expressed in terms of integrands with eight or less denominators. The 
recursive procedure described in Section 2.2 leads to the following multipole decomposition 



X, 



Tl 



=1 Dil "' Dis i, «ir=l Dil "' Dl7 



i 1 «i 8 =i - - - s n«i 7 =i " 11 _ *■' h<i 2 =i ~ 11 12 i=l 



where i a « % stands for a lexicographic order i a < i a+ \ < . . . < i^-i < %■ Equivalently, 
the numerator decomposition formula reads, 

n n n n 

Mi... n = yi a h-*8 n Dj + yi Ai i-^ n Dj+---+ 

il«ig=l jVii,...,i g ii«ij=l j^i 1: ...,i 7 

n nun n 

+ E II n > ■ Y, x \l l} ■ Q U n ■ ( 2 - 13 ) 

il<i2=l jVii,*2 *=1 jV^ i = l 

The residue A^...^ is obtained from the corresponding rank r^...^ integrand 2^...^ using 
the following procedure 

1. Decompose the loop momenta in using two bases {Tj}^!,...^ and {ej} J= i r .. ) 4, 



l 



A 1 



+ J>rf, fc" = -r£ + 5>e£. (2.14) 

i=i i=i 

In this case z = (j/i, . . . , 2/4, xi, . . . , X4). 

2. Consider a generic rank r^...^ polynomial in z 

M,.. ife ( Z )= £ ay(n*?j> ^) = {jeN 8 :^j 4 <r}. (2.15) 
JeJ(r h ... ik ) Vi=l / i=l 

3. Choose a monomial order and construct a Grobner basis Gi x ---i k = {51 (z), . . . , g m (z)}, 
generating the ideal Ji x ...{ k = {D h ,. . . , D ik ). 



4. Divide A/^...^ modulo Qi v ..i k holding the remainder Aj r 



The integrand decomposition ( 2.12 ) allows to express the amplitude in terms of Master 
Integrals (Mi's), associated to diagrams with 8-, 7-,. . . , 2- denominators. Depending on the 
powers of the integration momenta appearing in the numerator, the multivariate division 
may also generate the single-cut residues A,, and the quotients of the last divisions, Q . 
These two contributions generate spurious terms only but they are needed for the complete 
reconstruction of the integrand. The term Q is non cut-constructible: its determination 
requires to sample the numerator away from the solutions of the multiple-cuts. 

We expect that the integrand reduction formula could be extended to d dimensions, 
where additional degrees of freedom related to X q and enter [65]. 

1 A potential ambiguity may arise in topologies with 9 denominators two of which are degenerate. How- 
ever in this case the one-loop sub-topology contains at least 6 denominators yielding thus a system of 
equations with no solution. 



- 8 - 




Di 
D 2 
D 3 
Di 
Ds 
D 6 
Lh 



= (k + p 2 ) 
= (k-pi) 

2 



(9- 



Pi) 



(g - P4) 

(9 -p* -ps) 2 

(5 + ft +P2 +P3) 



(a) 




Di 

a, 

D 3 
D 4 
D s 
D 6 
D 7 
D 8 



= fc" 

= (fc + P2 ) 2 

= (k + q - pa 
2 

= 9 

= (<7 + £>3) 2 

= (<7-P4) 2 

= (9- 

= (g- 



P5) 



Pi -Ps) 
- k + p2 + P3) 



(&) 




D 2 
-Da 
Di 
Ds 
D 6 
D 7 
D s 



g 
(9 



Pi) 



= (9 + P4) 
= k 2 

= (fc + P5 ) 2 

= (k-pi) 2 

= (ft - g + p 3 +p 5 ) 2 

= (k-q-pi - Pa) 2 



(c) 



Figure 1: Five-point diagrams entering the amplitudes in A/" = 4 SYM and N = 8 SUGRA. They 
are the pentabox diagram (a) the crossed pentabox diagram (b) and the double pentagon diagram 
(c). For each diagram, the definition of the denominators is shown as well. 



3. Five-point Amplitudes in N = 4 SYM 

We apply the integrand reduction to five-point two-loop scattering amplitudes in J\f = 
4 Super Yang-Mills (SYM). The five-point amplitudes get contributions from the three 
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diagrams depicted in Fig. |l[ We introduce the notation X^ 4 '^ ( AA 4 '^) to denote the 
integrand (numerator) of the diagram in Fig. m (i) in J\f = 4 SYM. The integrand of 
diagrams in Fig. [I] are 

L>i---±J S 

4*3 («, *) = M p {q n k) , A/?? k) = 2 g • Ul + ft , (3.2 

Zft? (3, *) = M £* iq n k) , A/ft? (ff, fc) = 2 g • tij, + 2 fc • t* + & + & , (3.3 
U\ ■ ■ ■ JJs 

The vectors u^, u^, and U3 and the constants are defined as [52] 

u l = J (735124(^5 - P3) + 734125 (f*4 - P3) + 745123 O5 ~ P4) + 2 712345 (^2 ~ Pi)) > ( 3 - 4 
U 2 = \ (723145(^2 - P3) + 724135 (^2 - P4) + 734125 (f>3 ~ P4) + 2 715234 (Pi ~ P5)) , (3.5 
U 3 = 7 (7l2345(Pi - P2) + 725134 (f>2 - P5) + 715234 (Pi - P5) + 2 734125 (P3 ~ P4)) , (3.6 



A = 7 (735124 O.S! + -7l.2 + - s :s7 + - 7'! 1. 1 2:> *J2 + * ir.J2:S ('*:<! + 7! 2 + *:sr.) 



+2712345(^23 - S13) ) , (3.7 
P2 = -. ( ~ (723145 + 724135) S23 + 734125(«15 + S34 + 2 S23) - 2 715234 (S13 - S35)) , (3.8 



/?3 = ^((712345 - 725134) Sl2 + 715234(S34 + Sl5 + 2 S12) - 2734125(^13 - S14)) • (3.9 

where the kinematic invariants Sjj and the functions 7 read as follows 

Sij = (Pi+Pj) 2 = 2{Pi- Pj) (3.10 
. -( [12] [2 3] [3 4] [4 5] [51] \ 

712345 \ v [ 14 ][23](i2)(34) - [12][34](14}(23)y { h { ' 

In A/" = 4 SYM, given the simple form of the numerators, the multipole decomposition 
of the integrands only requires one iteration. The numerators can be decomposed as: 

7 

Nh-'^iV' k ) = A 12345678 + ^2 ^l-(i-l)(i+l)-8 A > X = a,b,C . (3.12) 

i=l 

The number of 7-ple residues of the integrands AA 4 ' a ) and AA 4 '^ is almost halved since the 
numerator depends on q only, thus A 1 ...(j_ 1 vj +1 v.. g = for z 7^ 4, 5, 6, 7. 

In the next subsections, we list the parametrization of the residues entering our com- 
putation, namely of the residues in Eqs. (|3.12j ). All the eightfold residues are related to 
maximum cuts. According to the maximum cut theorem [44], the number of coefficients 
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Figure 2: Cut-diagrams of the sevenfold cuts. Starting from the left, we show the diagram of the 
cut (1234568), (1234678), (1234578), and (1235678). 



cut 


bases 


z 


Monomials in the residue 


(12345678) 


Eq. Q3.13D 


(y4,y3,y2,yi,x4,x 3 ,x 2 ,x 1 ) 


Sl2345678 = i 1 > x l , 2/1 , y 2 } 


(1234568) 


Eq. (|3.13) 


{y4,V3,y2,yi,X4,X3,X 2 ,Xl) 


^1234568 = {l,Xl,x\,x\,x\,X 2 ,XxX2, 

xlx 2 ,x 3 1 x 2 ,y 1 ,x 1 y 1 ,xly 1 ,x 3 1 y 1 ,xfy 1 , 

2 3 2 2 

x 2 yi, x\x 2 y\, x x x 2 yx, x 1 x 2 yi,y 1 ,x\y l) 
X2yf, 2/i , x%yf, x 2 y\, yf, xiyf, x 2 y$, y 2 , 
xiy 2 , 2/12/2, 2/12/2,2/12/2} 


(1234678) 


Eq. (|3.13) 


{y4,V3,y2,yi,x4,x 3 ,x 2 ,xi) 


>5i234678 = ^1234568 


(1234578) 


Eq. (|3.14) 


{Vi,V3,y2,yuXi,X3,X2,Xl) 


<Sl234578 = {I, X 2 , X 2 , X 2 , X 2 , X4, X 2 X4, 

x 2 Xi, xf,x 4 , y 1 ,x 2 y 1 ,xly 1 ,xly 1 ,x%y 1 ,x 4 y 1 , 
x 2 x i y 1 ,xlx i y 1 ,x 3 ,x 4 y 1 ,yf,x 2 yf, x 4 yf, yf, x 2 y\, 
xavI, 2/i , x 2 y\, x 4 yf, 2/4, x 2 y 4 , y x y 4 , 

yiV4, 2/12/4} 


(1235678) 


Eq. Q3.15) 


{y4,y3,y2,yi,x4,x 3 ,x 2 ,x 1 ) 


^1235678 = {l,Xi,xl,xl,xf,X 2 , 

xix 2 ,xlx 2 ,xfx 2 , j/i, xxyi, xfyi, 

xlyi,xlyi,x 2 yx,xix 2 yi,xlx 2 yi,xlx 2 yi, 

yi, xiyf,x 2 yl,yf, xiy\, x 2 y\, 

yf, xiyf, x 2 y\, y 3 , xiy 3 , 2/12/3, 2/i2/3, 2/i2/3} 



Table 1: Set of monomials parametrizing the residues entering the decomposition of the five-point 
pentabox diagram. They have all been found using degree lexicographic monomial ordering. For 
each cut the bases and the chosen ordering for loop-variables are shown as well. 



needed to parametrize the residue of a maximum cut is finite and equal to the number 
of the solutions of the corresponding cut. We found the most general parametrization of 
the eightfold residue, which is process-independent and valid for numerators of any rank 
in both q and k. The parametrization of the sevenfold residues is also process- independent 
and it is given for the case of renormalizable numerators of rank six at most, which is 
more than we need for the applications presented in this paper. The parametrization of 
higher rank numerators can be obtained including additional terms with higher powers of 
the loop momenta in the residues [53] . The number of coefficients required to parametrize 
the residues agrees with an independent computation performed using a technique based 
on Gram determinants. 
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Figure 3: Cut-diagrams of the sevenfold cuts. Starting from the left, we show the diagram of the 
cut (1234568), (1234678), (1234578), and (1235678). 



3.1 Residue of the planar pentabox 

The decomposition of the pentabox diagram in Fig. |l| (a) requires the parametrization of 
the residues of the eightfold cut (12345678) and the sevenfold cuts depicted in Fig. ^. The 
parametrization is obtained using the procedure described in Section |2| The relevant bases 
are 



r = 


0", et = 


A< 

Ps> 


A 4 A* 

e 2 ~~ Pa 


j> <3|V|4] _„ (4| 7 *|3] 
3 — 2 ' ~~ 2 ' 


< K = 


M , rf = 


A 1 


a a 

t£ = Vx 


_ M <2| 7 "|1] _ M (1| 7 "|2] 
3 ~ 2 ' ' 4 2 ' 


Xl = 


(g-pi) 

[PVP2) ' 


™ r _ (9'Pa) 
X2 ~ (P1-P2) 


_ (fc-p 4 ) _ (fc-pa) . 

' yi (P3'P4) ' yZ (P3'P4) ' 


r = 




A 1 

Pi) 


e 2 — Pz 


(3|2|1](1| 7 ^|3]±(1|2|3](3| 7 ^|1] 

3,4 — | 




0", rf = 


A 1 
Pi' 


a fi 
T 2 = P3 


_ M (3|2|1](1| 7 ^|3]±(1|2|3](3| 7 ^|1] 
T3 ' 4 ~ 4 


Xl = 

V 


(<?'P3) 

(pi-Ps) ' 


X,i 




2/1 ~~ (P1-P3)' ^ 4 ' 


r = 






(3|2|l]<l|y|3]=F(l|2|3]<3|7"|l] ^ „ „ „ 

I e 2 — P3 ' e 3 — Pi 


A 1 

' Po = 


A 4 A 4 
— Pii r l,4 




(3|2|1](1| 7 ^|3]t(1|2|3](3| 7 ^|1] _ M ^ _ M ^ 

I J T 2 — ^3 ' T 3 — Pi 


Xl = 

V 


((g-p 4 )-ei) 






P4)'Pl) _ On) _ (fc-ps) 
l'P3) ' ^ t{ ' ^ 3 (P1-P3)' 



(3.13) 



(3.14) 



(3.15) 



The residue of each cut (ii . . . i K ) is written in terms of a set of monomials *Sj r 
in Table § 



collected 



3.2 Residue of the crossed pentabox 

The diagram in Fig. [I] (b) is decomposed in terms of the residue of the eightfold cut 
(12345678) and of the residue of the sevenfold cuts in Fig. |3[ Each residue can be expressed 
in terms of a set of monomials, as shown in Table The parametrization is obtained using 
the multivariate polynomial division described in Section |2[ 
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cut 


bases 


z 


Monomials in the residue 


(12345678) 


Eq. (|3.13) 


(^4,^3,^2,2/3,^4,^1,^2,2/0 


^12345678 = {1,^1,2/1,2/2} 


(1234568) 


Eq. (|3. 13) 


(y4,y3,y2,yi,x4,x 3 ,x 2 ,x 1 ) 


^1234568 = {l,aci,a;f,ij,xf,a;i,x;,a:2) 

^1^2, X±X2, xfx2, xfx2, x\x2, 

y 1 ,x 1 y 1 ,xly 1 ,xly 1 ,xfy 1 ,xly 1 ,x 2 y 1 , 
xxx 2 yi , x\x 2 yi , x\ x 2 yi , x\ x 2 yi , y\ , x x y\ , 
X2yl, 2/i , xiVi, x 2Vi, 2/i, xiyf, x 2 yf, 2/2, 
2:12/2,2/12/2,2/12/2,2/12/2} 


{LZoQO to} 


TTn ft! Til 


{y4,V3,y2,yi,X4,X3,X2,Xl) 


c c 

01234678 — '-'1234568 


(1234578) 


Eq. (p. 14) 


(y4,y3,y2,yi,x 4 ,x 3 ,x 2 ,x 1 ) 


<?1234578 = {1, x 1i X21 X 2 , X 2 , X 2 , X 2 , £4, 

x 2 Xi, x\x 4 , x\x 4 , x%X4, x\x i ,y 1 ,x 2 yi, 
x\y\ , x\y\ , x\y x ,x\y\, x 4 y x , x 2 x 4 yi , x\x 4 y x , 
x\x 4 y x , x\x 4 yi ,yl,x 2 yf, x 4 y\ ,yl,x 2 yf, 

x4Vi, 2/1, x 2 yf, x A yf, y 4 , X2V4, 2/12/4, 2/12/4, 2/12/4} 


^IZoOD / o J 


riq. (ji.lo) 


(2/4, 2/2, 2/3, 2/1, ^4, ^3, #2, ^l) 


c _ /It ^,2 3 4 5 6 

>- 5 1235678 — l 1 , #1, Xi, X x , X x , X 1; X % , X2, 

X\X2,X\X2,X\X2, x\x2,x\x2,y\ 1 

xiyi,xlyi,X2yi,xix 2 yi,yl,x 1 y1, x\y{, 
Z22/?, xix 2 y\, y\, xiyf, x\y\, x 2 y\, xix 2 yf, 
2/i, xiyf, xjyf, x 2 yf, xix 2 yt, 2/3, x 4 y 3 , 
2/12/3, 2/i2/3,2/i2/3} 



Table 2: Same as Table [j], but for the five-point crossed pentabox diagram. 



3.3 Residue of the double pentagon 



The decomposition of the double pentagon diagram requires the parametrization of the 
residues of the eightfold cut (12345678) and all the sevenfold cuts. However, the topology 
in Fig. [I] (c) is invariant under the transformation 



Pi <r> P3 , 



Pa Ps > 



^ o A;^ 



(3.16) 



thus the only sevenfold cut needed are (1345678), (1245678), (2345678), and (1234567), 
depicted in Fig. ||. The remaining sevenfold cuts can be obtained using the transforma- 
tion (|3. 16 ) . The eightfold cut is a maximum cut. It exhibits eight solutions and it is 
parametrized by eight coefficients, in accordance with the maxim-cut theorem. 



The sets of monomials parametrizing the relevant residues are collected in Table |3[ 
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cut 


bases 


z 


Monomials in the residue 


(12345678) 


Eq. (|3.17) 


(2/4,2/3,2/2,2/l,x4,x 3 ,x2,xl) 


^12345678 = {1, 2/l, Xl , 2/1^1, X?, X?, ^2, XiX 2 } 


(1345678) 


Eq. (|3.17) 


(x4,x 3 ,x 2 ,xi, 2/4, 2/3, 2/2,2/1) 


<$1345678 = {1,2/1,2/1,2/1,2/1,2/1,2/1,2/2,2/12/2, 

o 4 c 

2/i2/2,2/i2/2,2/i2/2,2/i2/2,a;i, 
2/iXi, 2/1X1, 2/1 £1, 2/iXi, 2/fxi, 
2/2^1, 2/12/2^1, 2/i 2/2X1, ViViXx, 
j/i 2/2X1 , x? , 2/1X1 , 1/2X1 , x? , 2/1 x? , 
2/2X?,xf,2/ixf,2/ 2 xf, 

X 2 , 1/1X2, XiX 2 ,XiX2,xfx 2 } 


(1245678) 


Eq. (|3.17) 


(x4,x 3 ,x 2 ,xi, 2/4, 2/3, 2/2,2/1) 


»$1245678 = ^1345678 


^OIOU 1 o J 


£jq. [p. La) 


(xi, X 3 , X 2 , X4, 7/3, J/4, y2,2/l) 


52345678 = {1, 2/1, 2/?, 2/?, 2/f, 2/l, 2/?, 2/4, 2/12/4, 

2/i2/4, 2/i 2/4, 2/i2/4, 2/l2/4, x 2 , 2/1X2, 

X4, 2/1X4, 2/1X4, 2/1X4, 2/1X4, 2/1X4, 

I/4X4, 2/12/4X4, 2/1^4X4, 2/i2/4^4, 

2/^2/4x4, x 2 x 4 , x|, 2/1 x|, 2/ 4 x|, 

X2X4, xi, 2/1X4*, 2/4X4, x 2 x^, x|, 2/1X4, 2/4xl| , 


(1234567) 


Eq. (p.17) 


(£4,2:3, £2, £1,2/4, 2/3, 2/2, 2/i) 


<$123567 = {1,2/1, 2/?, 2/1, 2/1, 2/2, 2/12/2, 2/l2/2, 2/i2/2, 

xi , 2/iXi , 2/iXi , 2/i^l , yfxi , 2/2X1 , 2/12/2X1 , 
2/i 2/2x1 , 2/12/2X1 , xf , yixf , 2/2X? , x? , 
j/i xf , 2/2X? , xf , 2/1 xf, 2/2xf , x 2 , 2/1X2, 

XiX 2 ,XiX2,XiX 2 } 



Table 3: Same as Table [|, but for the five-point double pentagon diagram in Fig. [I] (c). 



5 4 5 5 4 5 




Figure 4: Cut-diagrams of the sevenfold cuts. Starting from the left, we show the diagram of the 
cut (1234567), (1245678), (1345678), and (2345678). 



They are obtained by multivariate polynomial division using the following bases 



r£ = 0", e£ 



xi 



A 4 

Po 



,r 2 



P4, 



Ps> 



P3' 



Pl> 



(g-pi) 

(P5-Pl) 



-P3. r i 

((g+Pa)-Pi) 

(P1-P3) 



X2 



(£Ps) 
(P5'Pi) ' 



Pi. 


c 2 


A 1 

= Ps 


Pi, 


T 2 


A 1 

= P3 



yi 



'3,4 



'3,4 



(4^3] 
2 

(5|7 M |1] 
2 

(fc"P3) 



(3|7^|4] 
2 

(1| 7 A*|5] 



y2 



(fc-P4) . 
(P3-P4) ' 



(P3-P4) ' 

(1|4|3](3|7^|1]±(3|4|1](1| 7 ^|3] 
4 

(1|4|3](3| 7 ^|1]±(3|4|1](1| 7 ^|3] 



X4 



((g+P3)- r 4) 
• 



yi 



(fc-P3) 
(Pl'P3) ' 



2/1 



14 



(3.17) 



(3.18) 



4. Semi-numerical integrand reduction 



In the previous section we illustrated how to determine the general structure of the residues 
by means of the multivariate polynomial division. Knowing this structure, we can pro- 
ceed and numerically perform the integrand reduction to extract the values of all process- 
dependent coefficients which appear in the residues. The decomposition can be checked by 
verifying the identity between the original numerator and its reconstruction, i.e. between 
l.h.s. and r.h.s. of Eq. ( 3.12| ), for arbitrary values of the integration momenta q and k. 



This procedure is known as global N = N test of the integrand reduction. 
4.1 Planar pentabox diagram 

Eightfold cut. The residue of the eightfold cut can be parametrized using the monomials 
in Table 

^12345678 = C12345678, + C12345678, 1 (<7 ' Pi) + C12345678, 2\k • P3) + C12345678, 3 {k • Pi) • (4.1) 

The number of solutions equals the number of coefficients, in accordance with the maxi- 
mum cut theorem. Therefore the four coefficients appearing in Eq. ( |4.1| ) can be obtained 
by sampling the numerator on the four solutions of the eightfold cut, where the decompo- 



sition ( 3.12 ) becomes 

= ^12345678 • (4.2) 

In our case we find that only C12345678, and C12345678, 1 are non- vanishing. 

Sevenfold cut. The residue of the generic sevenfold cut (i\ ■ ■ ■ ij) appearing in Eq. (|3.12j ) 
can be parametrized using the results listed in Table |l[ For the process at hand, the 
structure of the numerator ensures that residue can be parametrized just by the a constant 
term, 

Ah-ir = <Hi-i 7 ,0 • (4-3) 

The actual value of Ci 1 ...j 7) o is obtained by sampling the numerator and the residue of the 
eightfold cut in correspondence of one solution of the sevenfold cut, where 

a _ M^il^k) - Ai2345678 

^ii-i 7 ~ 778 n • 

LLh^ii,...i 7 u h 

The multipole decomposition of the integrand I^!^ becomes 

T-(4,a)/ ,\ _ c 12345678,0 + ^12345678, l(g ' Pi) , c l-(t-l)(i+l)-8, , . -s 

m + h \\UD h ' <45) 

This result also shows the decomposition of the integral as linear combination of two Mi's 
with eight denominators and four Mi's with seven denominators. 
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4.2 Crossed pentabox diagram 

Eightfold cut. The residue of the eightfold cut is parametrized as (cfr. Table 0) 

^12345678 = C12345678, + C12345678, 1 (<? ' Pi) + C12345678, 2(k • P3) + C12345678, 3$ " Pa) ■ (4-6) 

The coefficients are obtained sampling the numerator at the four solutions of the maximum 
cut (12345678), where 

= A 1234 5678 • (4.7) 

The only non- vanishing coefficients are C12345678, and C12345678, i- 

Sevenfold cut. The generic sevenfold cut appearing in the multipole decomposition of 
lt ,b 8 is (k ■■■17) G {(1234568), (1234578), (1234678), (1235678)}. The structure of 
guarantees that the only non- vanishing coefficient is the one of the monomial 1, i.e. 

Aj r ..i 7 = Cj r ..j 7) o • (4.8) 

We sample the numerator and the residue of the eightfold cut at a solution of the cut 
(ii ■ ■ ■ ij) and we get c^.-.j^o by using the relation 

a _ A(£f-M0 - ^12345678 Q , 

LLh^ii,...ir u h 

Within the numerical precision, the coefficients of the crossed pentabox turn out to be 
equal to the ones of the planar pentabox. This is expected since the two diagrams share the 
same numerator and also the denominators appearing in its decomposition are in common 
between the two. 

The integrand xj 4 '^ is decomposed as follows: 



-7-(4,6)/ ,s _ C12345678.0 + C12345678, l(g ' Pi) , ^ c l-(t-l)(t+l)-8, (Ain\ 

W**>- o^w s + ^ x\UD h ■ (410) 

As for the previous diagram, this result yields the decomposition of the integral as linear 
combination of two Mi's with eight denominators and four Mi's with seven denominators. 

4.3 Double pentagon diagram 

Eightfold cut. The parametrization of the residue of the eightfold cut is given in Table || 
and can be written as 

^12345678 = c 12345678, + c 12345678), l(& ' Pz) + c 12345678, 2{Q • Ps) 

+ C12345678, 3{l ' Pi) + C12345678, a{Q ■ Psf + C12345678, h{l ' P5)(Q ' Pi) 

+ C12345678, &(k • Pz){q ' Ps) + C12345678, 7(<? ' Pbf ■ (4.11) 

The eightfold cut is a maximum cut, thus the eight solutions of the eightfold cut allow to 
determine the coefficients in Eq. ( |4.1l| ) using the relation 

J^S = A 123 45678 , (4.12) 

which holds at the solutions of the eightfold cut. The non- vanishing coefficients are 

C12345678, i for % < 4. 
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Sevenfold cut. The numerator has rank one and it is easy to see that the generic sev- 



enfold residue A^...^ entering Eq. (3.12) can be parametrized by a constant term, i.e. 

A h ... i7 = • (4.13) 

The value of c^...^^ is obtained by sampling the numerator and the residue of the eightfold 
cut at a solution of the cut (i\ ■ ■ ■ i?), where the relation 

a _ Afj 4 :.l\q,k) - Ai 234 5678 ,x 

llh^h,...ir u h 

holds. The multipole decomposition of the integrand of the double pentagon reads as 
follows 

r(4,c)/ M _ c 12345678,0 + C(i2345678, l(^ ' P3) + C12345678, 2(9 • P5) + C12345678, 3(9 ' Pi) 



D 1 ---D 8 

+ E Cl -t" ( T"° - (4 - 15) 

i=l llh& D h 

The corresponding decomposition of the integral is a linear combination of four Mi's with 
eight denominators and eight Mi's with seven denominators. 

5. Five-point amplitudes in A/ = 8 SUGRA 

In this section we decompose the five points two-loop amplitudes in M = 8 Super Gravity 
(SUGRA). They get contribution from the three diagrams in Fig. |j. In the following we 
denote the integrand (numerator) of the diagram in Fig. (t) by T&ti ( AA 8 '^). The 
integrands in M = 8 SUGRA are obtained by squaring the corresponding integrand in 
A" = 4 SYM [52]. 

The numerators are of rank two in the loop momenta. Following the same machinery 
as in the case of N = 4 SYM, we show that their decomposition can be expressed in terms 
of 8-, 7-, and 6-denominator integrands, 



Ml 8 : X 8 \q,k) = A 123 45678 + E 



i=l 

8 

+ A i-(i-i)(i+i)-(i-i)0'+i)-8A-Dj , x = a,b,c. (5.1) 
i<j=i 

The corresponding decomposition for the integrands xj 4 '^ , xj 4 '^ and zj 4 '.^ reads 

Ui,x), , N Ai2345678 . A l-(i-l)(j+l)---8 

+ Y~ Ai-.( i -i)( i +i)...b-i)0 + i>.- g _ x = aAc , (5. 2) 

i<j=i n^ij D h 

Since the numerators A/"| 8 !g^ and M^'^ are of rank 2 in q and independent of k, their 
decomposition is significantly simplified. Indeed in these cases Ax...u-i)u+i)...s = for 
i ^ 4,5,6,7 and A!...^^...^^.^..^ = for i, j ^ 4,5,6,7. 
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5.1 Semi-numerical computation 

In this section we briefly describe the numerical decomposition of the three numerators. We 
checked the decomposition by verifying that at arbitrary values of q and k the numerator 
and its reconstruction are equal. 

Planar pentabox 

Eightfold cut. The computation of the residue of the eightfold cut, A12345678 follows 
the same pattern as the M = 4 SYM planar pentabox, described in Section 4.1. The 
non- vanishing coefficients are 012345678,0 and C12345678 i- 



Sevenfold cut. The residue of the generic sevenfold cut (if ■■17) in Eq. ( |5.1| ) can 
be parametrized in terms of the monomials collected in Table [I]. The structure of the 
numerator guarantees that the residue contains rank one terms at most, 

A-h-ir = Qi-i7,o + Cii-i7,i(<? + w o) ■ w i + <Hi-i 7 ,2(q + w ) ■ w 2 

+c il ... i7)3 (k + u> 3 ) • Wi + c ix ... ilA (k + w 3 ) ■ w 5 . (5.3) 

The momenta w^ depends on the cut (i\ ■ ■ ■ if). The actual value of the coefficients can be 
obtained by sampling on five independent solutions of the sevenfold cut, where 

a _ N^'.^iq, k) - Ai2345678 ^ 

LLh^ii,...ir U h 

In this case, all the coefficients but Cj 1 ...j 7) o, ^1-47,1, and cn—i 7 ,2 vanish. 

Sixfold cut. The numerator has rank 2 in q and is independent of k. The 

only non- vanishing term in the residue of the generic sixfold cut (i\ - ■ -i§) in (5.1) is the 
constant, therefore we have 

^ii—ie = Cii— ie,0- (5-5) 

The actual value of the constant can be obtained evaluating the decomposition (^^) at 
one solution of the sixfold cut, where 

a _ ■^l. 8 .8^(g' fc ) ~ ^12345678 -V •• /.• ■/,, -v 

llh^i 1 ,...i 6 1J h h=4 u h 

h^i 1 ,...i (j 

After polynomial fitting of A 12345678 > \i—i 7 and Aj r ..j 6 , the resulting multipole de- 
composition of Eq. (5^) contains 20 non- vanishing coefficients two of which are spurious 
(i.e. their contribution vanish upon integration) while the others give rise to Mi's, namely 
two with eight denominators, ten with seven denominators and six with six denominators. 

Crossed pentabox 

As in the M = 4 SYM case, the crossed pentabox in Fig |l] (b) has the same numerator and 
the same decomposition as the planar pentabox. Therefore the coefficients of the former 
are exactly the same as the coefficients of the latter. 
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Double pentagon 



Eightfold cut. The computation residue of the eightfold cut of the double pentagon fol- 
lows the same lines of the M = 4 SYM double pentagon (see Section 4.3). The parametriza- 



tion of the residue is given in Eq. (4.11). In this case the only vanishing coefficient is 

C12345678, 7- 

Sevenfold cut. The residue Aj 1 ...j 7 of the generic cut (i\ •••£7) can be parametrized 



using Eq. ( |5.3|) . At the sevenfold cut - ■ ■ 17) the decomposition (pjj) reduces to 



■A/ffl jj^jQ ~ A 1234 5678 
Ylh^h,...ir Dh 



a _ ^ v l-8 Vj> "V ~ ^12345678 7 n 



The coefficients are then computed by sampling Eq. (5.7) at five solutions of the sevenfold 
cut. The non- vanishing ones are those multiplying constant or linear terms in the loop 
momenta. 

Sixfold cut. The residue A^...^ of the generic sixfold cut (if ■ ■ i§) can be parametrized 
by a constant, as in Eq. ( |5.5| ). The constant is computed using one solution of the sixfold 
cut and the expression of the decomposition ( |5.1[ ) at the sixfold cut, 

a ■Afi 8 .'.l > {q,k) - Ai2345678 -V • •••/••••<„ , , 

iv " i6 ~ ffl n ^ — Dl — ' ( } 

We find that A123456 = 0, while the residues of all the other sixfold cuts are non-vanishing. 

After polynomial fitting of Ai2345678> Aj r ..j 7 and Aj r ..j 6 , the resulting multipole de- 
composition of Eq. (5.2) contains in this case 74 non-vanishing coefficients, four of which 



are spurious. The integral can be decomposed as a linear combination of seven Mi's with 
eight denominators, 36 Mi's with seven denominators and 27 Mi's with six denominators. 



6. Analytic integrand reduction 

In this section we perform the reduction of the five-point amplitudes analytically. We apply 
a two- loop generalization of the integrand reduction through Laurent expansion formulated 
in [53]. As in the one loop case, the Laurent expansion allows to find simpler formulas 
for the coefficients entering the decomposition. Moreover, the subtraction of the higher 
residues can be performed at the coefficient level rather than at the integrand level. Indeed 
the Laurent expansion makes each function entering the reduction separately polynomial. 
Therefore the subtraction can be omitted during the reduction and accounted for correcting 
the reconstructed coefficients. For simplicity we will focus on the rank-one numerators in 
the five-point amplitudes of N = 4 SYM, however the method can be extended to higher 



rank numerators as described in Section 6A for the planar pentabox diagram in AT 
SUGRA. 
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6.1 Planar pentabox diagram 

The analytic decomposition of the pentabox diagram in Fig. [l] (a) we are about to discuss, 
is valid for any numerator of the type ( |3.1| ), i.e. for any rank-1 numerator depending on q 
only. Indeed our computation is carried out for generic u 1 ^ and f3\; the results for M = 4 
SYM will be recovered at the very end, using Eqs. (JO|) and (13.1 



Eightfold cut. The four solutions of the eightfold cut (12345678) are 

f {45) (3|7M|4] ( 32 ) (117 M 12] \ r „ ^ _ / (4 5) (3| 7 "|4] (15) (2| 7 "|1] 

l9l ' lj "V(35) 2 '(13) 2 J' i<?2,fc2j " V(35) 2 '(25) 2 

fjjg] <4|-y^|3] [15] i n r /[45] (4|Y*|3] [32] (2| 7 "|1] 

W3,3j "V[35] 2 '[2 5] 2 J' 194 ' * 4 J " V [3 5] 2 '[13] 2 



The general parametrization of the residue A12345678 is given in Eq. (|4.l| ). The simple form 
of the numerator Afit'.^ implies that the coefficients C12345678, 2 and C12345678, 3 vanish. The 
non- vanishing coefficients are obtained by sampling at the solutions (q%,ki) and (53, ^3) 
only. The outcome is 

1 

"(54)(31)[53][41]-(53)(41)[54][31] 

x ((5 4) (4 1) (3|«i |4] [5 4] [3 1] - (5 4) (3 1) (4|«i |3] [5 4] [4 1] 
-A (5 4) (3 1) [5 3] [4 1] + /3i (5 3) (4 1) [5 4] [3 1]) (6.1) 



C12345678, 



(5 4)(3|n 1 |4][5 3]-(5 3)(4|n 1 |3][5 4] 

C12345678, 1 ^ ^ ^ ^ ^ ^ ^ ^ _ ^ ^ ^ ^ ^ ^ ^ ^ ■ \y- ) 

Sevenfold cuts. We discuss the generic sevenfold cut (i\ ■ ■ ■ it) appearing the decompo- 
sition (3.12). For later convenience, we define the uncut propagator 

Di 8 (q,k) = (q + P ig ) 2 , with i 8 € {1, . . . , 8} and i 8 ^ h, ■ ■ ■ h- (6.3) 

The momentum Pf" is a linear combination of external momenta and it can be inferred 
from Fig. |l] (a), for instance Pj = —p^ — pg. The simplicity of the numerator allows to 
parametrize the residue using the constant term c^.-.j^Q, cfr. Eq. (|4.3|) . 

Every sevenfold cut of this topology exhibits a t-dependent family of solutions of the 
type: 

(«?>*?)= (<,!* + <,0><) • ( 6 - 4 ) 



The coefficient c^.^o can be computed evaluating Eq. ( [4.4|) at the solutions (|6 . 4|) . The 
computation can be simplified performing a Laurent expansion around t = oo, 



■A/l.^OZlA) Ai2345678(gl,fcl) 



— Cij -17, 

J t— >oo 



Di a (qi,h) D^q^ki 
Indeed in general neither 

12345678 (<?1, & 



(6.5) 



nor — — — (6.6) 



D i8 (qi,h) D is (qx,ki) 
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are polynomial in t but only their difference is. However, their truncated Laurent expansion 
obtained neglecting 0(l/t) terms is polynomial in t, namely in this case a constant 



^12345678 



A 8 (qi,h) 



= 6i 1 ...i 7 , + O( - 



(6.7) 



Therefore, similarly to the one-loop case [53], the coefficients 7ii 1 ...j 7 o an d ^i...i 7 ,o can be 
computed separately obtaining the coefficient c^.-.j^o by their difference, 



'ix— if, 



(6.8) 



The subtraction can be performed at the coefficient level rather than at the integrand level. 
Moreover the known structure of A12345678 allow to compute the coefficient &i r ..i 7l o once 
and for all, irrespective of the actual form of the numerator: 

c 12345678, ' v qi,l) 



hi-i r ,0 

The coefficients r%...i 7i o read as follows 

"1234568, " 



2(Pi 8 +v qu0 ) -v qul 



(6.9) 



"1234678.0 



"1234578, 



"1235678,0 



Ul 


V QlA 


P5 


%,1 


Ul 


V qiA 


P3 


V qiA 


Ui 


V qiA 


Pi 


V qi,l 


Ul 


V qi,l 


P4 


V qi,l 



with 
with 
with 
with 



J qi,l 



J qi,l 



J qi,l 



(3|7 M |4] 
2 

(5| 7 "|4] 



(317^345] 



(5|7 M 1^543] 



(6.10) 



Sbc 



-pZ 



(6.11) 



The massless momentum P a b c is defined as 

As already stated, the foregoing discussion applies to any numerator of the form given 
in Eq. ([0|). By using the explicit expressions of u\ and fi\ given in Eq.s ( p.4| ) and (|3"77|), 
we can write down the results for the coefficients in W = 4 SYM in terms of the functions 
7 defined in Eq ( [3.11 ) 

^(2712345(523 - si3 - S45) + 2734125 S12 

+ (735124 + 745123)(«12 + S34 + S35 + S45) 

— 2 712345 
1 



c 12345678, 



C12345678, 1 
C1234568, 

C1234578, 

C1234678, 

C1235678, 



("735124 — 745123 + 2 712345) 
("734125 + 735124 + 2 745123) 
(—735124 — 734125 — 2 712345) 
(-745123 + 735124 + 2734125)- 



(6.12) 
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The complete integrand decomposition is obtained plugging the coefficients of Eq. ( |6.12j ) 
in Eq. (|4.5| ). These results are in agreement with the ones found in the numerical compu- 
tation of Section 4.1. 



6.2 Crossed pentabox diagram 



As already noticed in Section |4.2| , the crossed pentabox of Fig. [l] (6) and the planar pentabox 
have the same decomposition. Indeed the numerator A/^g and the numerator A^ 4 ^, 



Eq. ( |3.2| ), are equal and are decomposed in terms of the same denominators, cfr. Eqs. ( 4.1 C| ) 



and (4.5) multiplied by ...Di 8 . Therefore the coefficients appearing the multi-pole 



decomposition ( 4.1C| ) are equal to the corresponding ones appearing the decomposition (| 



of the planar pentabox. 



6.3 Double pentagon diagram 

The numerator A^g , eq. (|3.3| ), is highly symmetric in q and k. Indeed the dependence 
on q and k can be disentangled and A^. 4 ]^ can be cast in the following form: 

M^(<i,k) = n^%)-T 1 {4 :l\k), (6.i3) 

where 

n^(t) = 2i ■ u 2 + ft, riti\k) = n?:i(k) \ P1 « P3 , P4 „ P5 (e.u) 



The vector u 2 and the constant ft are defined in Eqs. ( |3.6|) and (|3.9[) , respectively. Although 
in principle there are no additional issues in performing the reduction of the full numerator, 
the computation can be simplified by performing the reduction of 1Z^'. C g(q) only. The 
reduction of T^'g^q), and thus of the full numerator M^.'.g (q,k), is then obtained by 



means of the substitutions (3.16). The numerator TZ^'°g(q) depends on q only and can be 
decomposed as follows: 

3 

Kh'iil) = A 12345678 + ^2 A l-(i-l)(i+l)-8 A • (6-15) 

i=l 

Eightfold cut. The solutions of the eightfold cut (12345678) are eight and they can be 
used to compute the eight coefficient parametrizing the residue 

^12345678 = C12345678, + C12345678, 1 (k • P3J + C12345678, 2{l ' Pb) 

+ C12345678, 3(9 ' Pi) + C12345678, i{<! ' Pbf + C12345678, 5(<7 ' Ps)(Q ' Pi) 

+ C12345678, &(k ' Pz){<l ' P5) + C12345678, 7(<7 ' Pbf ■ (6.16) 

The rank of T^.'.'g implies that for i > 4 C12345678, i = 0. The simplicity of the numerator 
simplifies the computation even further. Indeed decomposing in the basis {pi}i=i,--- ,4, 
and using the conditions D% = D 2 = D3 = we get 

Cl2345678,0 = 02, C12345678, 1 = . (6-17) 
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The two remaining coefficients can be obtained by sampling the numerator on two solutions 
of the eightfold cut, e.g. 



(14)[41](45) + (3 5) [3 5] (4 5) - (4 2) (3 5) [2 3] (5[ 7 M |1] (3 5) (4| 7 "|3] 

(2 5) [12] (4 5) 2 ' (45) 2 

(23)[42][35] + (41)[41][45]-(35)[35][45] (1|t"|5] [3 5] (3| 7 ^|4] 



(21) [2 5] [4 5] I 

The missing coefficients read as follows: 

,( 4|l|3](3|n 2 |4]-(3ll|4](4|^|3] 

(4|5|3](3|1|4]-(3|5|4](4[1[3] ' 
(4|5|3](3|u 2 |4]-(3|5|4](4|tt 2 |3] 
(4|5|3](3|1|4]-(3|5|4](4|1|3] " 



' [4 5] 



C12345678, 2 — ~ 2 
C12345678, 3=2 



(6.18) 



Sevenfold cuts. We consider the generic sevenfold cut (i\ ■ ■ ■ ij) appearing the decompo- 
sition (6.15). The solutions of the cut can be cast into one-parameter families. In particular 
each cut allow for a solution with the following asymptotic behavior 



£ 1 / + lo + o(|),< 1 ,iH< )0 + oQ) fort^oo. (6.19) 



We compute the coefficient C1345678, by evaluating the decomposition ( p. 15 ) at one solution 
of the sevenfold cut, and by expanding around t = 00: 



12345678 



D i8 (qi,h) D is (qx,kx) 



"ii—ir, 



(6.20) 



t— >oo 



The denominator Di s is written in terms of P£ as in Eq. (|6.3| ). The actual form of Pj^ is 
inferred from Fig. [j] (c). Also in this case the t — >■ 00 limit makes both 

A 12345678 (<7l, h) 



Di s (qi,h) 



and 



Di 8 (qi,h) 



(6.21) 



polynomial in t, 



Di s (qi,ki) 



Ai2345678(9l)^l) 



,0 + O 



.22) 



The coefficients h^...^^ and 6i 1 ...i 7) o can be computed separately, and 

C»i— »7,0 = ^ii---i 7 , — &ii-i 7 ,0 • (6.23) 

Therefore the subtraction can be performed at the coefficient level via the universal function 

C12345678, 2(^91,1 " Pb) + C12345678, 3(^,1 ■ Pi) 



bii—iy, ( 



2(Vg lfi +Pi 8 ) 



(6.24) 
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The coefficients h^...^^ are given by 



^1345678, 



^1245678, 



^1234578, 

We define: 



u 2 ■ 






Vqi 


_ u 2 ■ 


% 


PA ■ 


% 


_ u 2 




P3 





m 
m 

113 

m 
m 



with Vq\ 



with 



with 



(5 2) [41] 
(4 2) [3 4] ' 
(5 2) [31] 
(4 2) [3 4] ' 
(3 5) [12] 
(3 4) [3 2]' 
(4 5) [12] 
(3 4) [3 2] ' 



(4|7 M |3] „ 

»7l n 1" 12?4, 



m- 



2 

(4| 7 "|3] 



W3- 



<3|7^|4] , <4| 7 "|3] 
o H??7 ; 



-. (6.25) 



-0"! + Vgl ~ 4(720-3 
2ct 2 

2(T2(Tg — CTl(T4 + 0"4-y/ CT^ — 4o"20"3 
2(T2(T5 



1o 2 07 — CTi(Tq + 0§\J a\ — 4o"2CJ3 
2CT2(T5 



(6.26) 



in terms of 



Oi = - 4(<7 7 cr 4 + 0- 6 cr 8 )(pi • p 2 ) , 
cr 2 = - 4(cr 6 cr 4 - of)(pi ■ p 2 ) , 
£J 3 = - 4o- 7 CJ 8 (pi -p 2 ) , 

0-4 = 2(pi • p 5 )(4 5) [3 1] - 2(p 2 • P 5 ) (4 5) [3 1] + (3 5) (4 2) [3 1] [3 2] - (4 5) (4 2) [3 2] [4 1] , 
0-5= -(3 5} (4 2) [3 2] [41] + (3 2)(4 5}[3 1] [4 2] , 

o-e = 2(pa, -ps) (3 5) [4 1] - 2(p 2 • p 5 ) (3 5) [4 1] + (3 5) (3 2) [3 1] [4 2] - (3 2) (4 5) [4 1] [4 2] , 
a 7 = (3 5) (5 2) [41] [12], 

erg = (4 5) (5 2} [31] [12]. (6.27) 
This completes the reduction of a generic numerator T^jt'.g of the form given in 



Eq. ( |6.14 ). As previously stated, the reduction of the full numerator JV^g of M = 4 



SYM can be recovered from the one discussed in this section, by means of Eq.s flglg ) 
and flP3) . To that purpose, we observe that the substitutions (|3,16j ) give the one-to-one 
mapping between denominators 

Dx o £> 4 , L> 2 <-> L> 6 , £> 3 O As, £>7 <-> Aj- 
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Putting everything together and using the definitions of u 2 and of Eqs. (|3.6| ) and ( 
we recover the full multipole decomposition of Eq. (4.15), with coefficients 



C1345678, 
C1245678, 
C1235678, 
C1234678, 
C1234578, 



Cl2345678,0 = J (^734125 (2«13 ~ 2«35 + 2s 2 3 + Sl5 + S34) 
~ 715234 (2S13 - 2s 35 + 2S12 + S15 + S34) 

- (723145 + 724135) S23 ~~ (725134 - 712345) («35 - «14 - Sl2)J 
C12345678, 1 = —2 734125 

C12345678, 2 = ^(2734125 _ 723145 _ 724135 _ 2715234 + 725134 ~ 712345) 
Cl2345678,3 = -(^734125 _ 723145 _ 724135 + 2715234 + 725134 ~ 712345) 
C2345678.0 = ^(2734125 ~ 723145 + 724135) 



( — 3734125 + 2723145 + 724135 _ 725134 + 712345) 
(734125 — 723145 — 2724135 + 725134 ~ 712345) 
( — 2715234 — 725134 — 712345) 
( — 2734125 + 715234 + 725134) 
(2734125 + 715234 + 712345) 



Cl234568,0 — ^(2734125 + 725134 ~ 712345) 
1 

C1234567, — ^( — 2734125 — 725134 + 712345)- 



(6.: 



The coefficients of Eq. fl6.28| ) enter the integrand decomposition of Eq. ( [4.15 ). These 
results are in agreement with the numerical computation of Section |4.3|. 



6.4 Higher-rank integrands 

The analytic reduction via Laurent expansion can be extended to numerators of higher 
rank. As an example we perform the computation of the sevenfold residues for the M = 8 
SUGRA two-loop five points planar pentabox, whose numerator has rank two in the loop 
momentum q. 

The most general parametrization of the residue of the generic sevenfold cut (if- if) 

is 



A il ... ir (q, k) = Ci v ..i 7t o + <v.. i7i i(g + w ) ■ w 1 + Ci v .. iT ^(q + w ) ■ w 2 . 



(6.29) 



The momenta Wq, w^, and w 1 ^ are a linear combination of the external momenta and 
depend on the cut (if ■ ■ ij). Their actual form is not relevant for this discussion but can 
be inferred from Table We consider two i-dependent solutions of the sevenfold cut: 



1i 1 



<,i* + <,o, < with =1,2 



(6.30) 
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The coefficients in Eq. (&2§j) can be computed evaluating ( |5.1|) at the solutions ( |6.30[) 

Ni S .'.8 (Qh k) ~ Ai2345678(%,^) _ « 



(ft) 



^■il'-ijVlii k) 
Cii-i r ,0 + Cii-ir, l(%,0 + ^o) • ^1 

+Cji-i 7 ,2(%,0 + ^o) • ^2 + Ci 1 ... i7j l(%,l • 

+Ci 1 -i T ,2(v qijl ■ w 2 )t, (6.31) 



where Dj g is defined in Eq. (|6,3j). The Laurent expansion around t = oo simplifies the 
computation. Indeed in this limit both 



and 



A 



12345678 



(ft) fej 



have the same polynomial structure of the residue, 



(6.32) 



•A/£?(«,*i) 



-^ig (ft) &i) 
A 12345678 (ft) 



■^ts (ft ! 



t— >oo 



t— >oo 



%-ir,0 + u 



(6.33) 



The expression of the coefficients is obtained by plugging the expansions ( 6.33| ) in Eq. ( 6.31| ) 
and by comparing both sides. In particular c^...^^ and c^...^^ are the solution of the 
system 



c h ...i 7A (v qi • wi) + c il ..^ 7>2 (% • ufe) = n^Li 7; i 

Cii-ir,i(% ' wi) + c il ...i 7 ,a(« ga ■ u> 2 ) = n !?-i 7 ,i 
The coefficient ^...j^o is given by 

Cii-ir^ = "<iii..i 7 ,o - V-i 7 ,0 ~ c 'r-<7,iKi,0 + ^o) ■ Wl ~ Qi-ir, 2 Ki,0 + U>o) • ^2 



(6.34) 



n [2] _ ^[2] 

in terms of the functions 



fti-i 7 ,l(%i,0 + w ) ■ Wi - Ci x ... i 7t 2(v q2 n + wo) • w 2 , (6.35) 



Cl2345678,l(Pl " % 



ii—ir, 



(6.36) 



[i] 



2 (%,o + -Pi 8 ) 

Eq. (6.35) shows that the coefficient c^...^ o can be written as the constant term 



n i —if o °^ the Laurent expansion of the integrand, corrected by two kinds of contributions. 
The first, , implements the the eightfold-cut subtraction as a correction at the 

coefficient level. The other terms are proportional to the higher rank coefficients of the 



same cut found as solutions of the system in Eq. (6.34). 



- 26 - 



7. Conclusions 

We recently proposed a new approach for the reduction of scattering amplitudes [44], 
based on multivariate polynomial division. This technique yields the complete integrand 
decomposition for arbitrary amplitudes, regardless of the number of loops. In particular it 
allows for the determination of (the polynomial form of) the residue at any multi-particle 
cut, whose knowledge is a mandatory prerequisite for applying the integrand- reduction 
procedure. We have also shown how the shape of the residues is uniquely determined by 
the on-shell conditions and, by using the division modulo Grobner basis, we have derived 
a simple integrand recurrence relation generating the multi-particle pole decomposition for 
arbitrary multi-loop amplitudes. 

In the present paper, we applied the new reduction algorithm to the two-loop five-point 
planar and non-planar scattering amplitudes in J\f = 4 SYM and M = 8 SUGRA in four 
dimensions, whose numerator functions contain up to rank- two terms in the integration- 
momenta. We determined all polynomial residues parametrizing the cuts of the correspond- 
ing topologies and sub-topologies. At the same time, the polynomial form of the residue 
defines the integral basis for the amplitude decomposition. For the considered cases, we 
identified independent integrals with eight, seven, and six denominators. 

Our presented approach is well suited for a semi-numerical implementation. The math- 
ematical framework it is based on is very general and provides an effective algorithm for 
the generalization of the integrand reduction method to all orders in perturbation theory. 
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